Transcriptomic analysis of HEK293A cells with a CRISPR/Cas9-mediated TDP1 knockout.

Tyrosyl-DNA phosphodiesterase 1 (TDP1) is a human DNA repair protein. It is a member of the phospholipase D family based on structural similarity. TDP1 is a key enzyme of the repair of stalled topoisomerase 1 (TOP1)-DNA complexes. Previously, with the CRISPR/Cas9 method, we obtained HEK293A cells with a homozygous knockout of the TDP1 gene and used the TDP1 knockout cells as a cellular model for studying mechanisms of action of an anticancer therapy. In the present work, we hypothesized that the TDP1 knockout would alter the expression of DNA repair-related genes. By transcriptomic analysis, we investigated for the first time the effect of the TDP1 gene knockout on genes' expression changes in the human HEK293A cell line. We obtained original data implying a role of TDP1 in other processes besides the repair of the DNA-TOP1 complex. Differentially expressed gene analysis revealed that TDP1 may participate in cell adhesion and communication, spermatogenesis, mitochondrial function, neurodegeneration, a cytokine response, and the MAPK signaling pathway.

A Knockout of Poly(ADP-Ribose) Polymerase 1 in a Human Cell Line: An Influence on Base Excision Repair Reactions in Cellular Extracts.

Base excision repair (BER) is the predominant pathway for the removal of most forms of hydrolytic, oxidative, and alkylative DNA lesions. The precise functioning of BER is achieved via the regulation of each step by regulatory/accessory proteins, with the most important of them being poly(ADP-ribose) polymerase 1 (PARP1). PARP1's regulatory functions extend to many cellular processes including the regulation of mRNA stability and decay. PARP1 can therefore affect BER both at the level of BER proteins and at the level of their mRNAs. Systematic data on how the PARP1 content affects the activities of key BER proteins and the levels of their mRNAs in human cells are extremely limited. In this study, a CRISPR/Cas9-based technique was used to knock out the PARP1 gene in the human HEK 293FT line. The obtained cell clones with the putative PARP1 deletion were characterized by several approaches including PCR analysis of deletions in genomic DNA, Sanger sequencing of genomic DNA, quantitative PCR analysis of PARP1 mRNA, Western blot analysis of whole-cell-extract (WCE) proteins with anti-PARP1 antibodies, and PAR synthesis in WCEs. A quantitative PCR analysis of mRNAs coding for BER-related proteins-PARP2, uracil DNA glycosylase 2, apurinic/apyrimidinic endonuclease 1, DNA polymerase ß, DNA ligase III, and XRCC1-did not reveal a notable influence of the PARP1 knockout. The corresponding WCE catalytic activities evaluated in parallel did not differ significantly between the mutant and parental cell lines. No noticeable effect of poly(ADP-ribose) synthesis on the activity of the above WCE enzymes was revealed either.

Usnic Acid Derivatives Inhibit DNA Repair Enzymes Tyrosyl-DNA Phosphodiesterases 1 and 2 and Act as Potential Anticancer Agents.

Tyrosyl-DNA phosphodiesterase 1 and 2 (Tdp1 and Tdp2) are DNA repair enzymes that repair DNA damage caused by various agents, including anticancer drugs. Thus, these enzymes resist anticancer therapy and could be the reason for resistance to such widely used drugs such as topotecan and etoposide. In the present work, we found compounds capable of inhibiting both enzymes among derivatives of (-)-usnic acid. Both (+)- and (-)-enantiomers of compounds act equally effectively against Tdp1 with IC50 values in the range of 0.02-0.2 µM; only (-)-enantiomers inhibited Tdp2 with IC50 values in the range of 6-9 µM. Surprisingly, the compounds protect HEK293FT wild type cells from the cytotoxic effect of etoposide (CC50 3.0-3.9 µM in the presence of compounds and 2.4 µM the presence of DMSO) but potentiate it against Tdp2 knockout cells (CC50 1.2-1.6 µM in the presence of compounds against 2.3 µM in the presence of DMSO). We assume that the sensitizing effect of the compounds in the absence of Tdp2 is associated with the effective inhibition of Tdp1, which could take over the functions of Tdp2.

Transcriptomic Analysis of CRISPR/Cas9-Mediated PARP1-Knockout Cells under the Influence of Topotecan and TDP1 Inhibitor.

Topoisomerase 1 (TOP1) is an enzyme that regulates DNA topology and is essential for replication, recombination, and other processes. The normal TOP1 catalytic cycle involves the formation of a short-lived covalent complex with the 3' end of DNA (TOP1 cleavage complex, TOP1cc), which can be stabilized, resulting in cell death. This fact substantiates the effectiveness of anticancer drugs-TOP1 poisons, such as topotecan, that block the relegation of DNA and fix TOP1cc. Tyrosyl-DNA phosphodiesterase 1 (TDP1) is able to eliminate TOP1cc. Thus, TDP1 interferes with the action of topotecan. Poly(ADP-ribose) polymerase 1 (PARP1) is a key regulator of many processes in the cell, such as maintaining the integrity of the genome, regulation of the cell cycle, cell death, and others. PARP1 also controls the repair of TOP1cc. We performed a transcriptomic analysis of wild type and PARP1 knockout HEK293A cells treated with topotecan and TDP1 inhibitor OL9-119 alone and in combination. The largest number of differentially expressed genes (DEGs, about 4000 both up- and down-regulated genes) was found in knockout cells. Topotecan and OL9-119 treatment elicited significantly fewer DEGs in WT cells and negligible DEGs in PARP1-KO cells. A significant part of the changes caused by PARP1-KO affected the synthesis and processing of proteins. Differences under the action of treatment with TOP1 or TDP1 inhibitors alone were found in the signaling pathways for the development of cancer, DNA repair, and the proteasome. The drug combination resulted in DEGs in the ribosome, proteasome, spliceosome, and oxidative phosphorylation pathways.

Biochemical Characteristics of iPSC-Derived Dopaminergic Neurons from N370S GBA Variant Carriers with and without Parkinson's Disease.

GBA variants increase the risk of Parkinson's disease (PD) by 10 times. The GBA gene encodes the lysosomal enzyme glucocerebrosidase (GCase). The p.N370S substitution causes a violation of the enzyme conformation, which affects its stability in the cell. We studied the biochemical characteristics of dopaminergic (DA) neurons generated from induced pluripotent stem cells (iPSCs) from a PD patient with the GBA p.N370S mutation (GBA-PD), an asymptomatic GBA p.N370S carrier (GBA-carrier), and two healthy donors (control). Using liquid chromatography with tandem mass spectrometry (LC-MS/MS), we measured the activity of six lysosomal enzymes (GCase, galactocerebrosidase (GALC), alpha-glucosidase (GAA), alpha-galactosidase (GLA), sphingomyelinase (ASM), and alpha-iduronidase (IDUA)) in iPSC-derived DA neurons from the GBA-PD and GBA-carrier. DA neurons from the GBA mutation carrier demonstrated decreased GCase activity compared to the control. The decrease was not associated with any changes in GBA expression levels in DA neurons. GCase activity was more markedly decreased in the DA neurons of GBA-PD patient compared to the GBA-carrier. The amount of GCase protein was decreased only in GBA-PD neurons. Additionally, alterations in the activity of the other lysosomal enzymes (GLA and IDUA) were found in GBA-PD neurons compared to GBA-carrier and control neurons. Further study of the molecular differences between the GBA-PD and the GBA-carrier is essential to investigate whether genetic factors or external conditions are the causes of the penetrance of the p.N370S GBA variant.

IPSC-Derived Astrocytes Contribute to In Vitro Modeling of Parkinson's Disease Caused by the GBA1 N370S Mutation.

Parkinson's disease (PD) is a neurodegenerative disorder that ranks second in prevalence after Alzheimer's disease. The number of PD diagnoses increases annually. Nevertheless, modern PD treatments merely mitigate symptoms rather than preventing neurodegeneration progression. The creation of an appropriate model to thoroughly study the mechanisms of PD pathogenesis remains a current challenge in biomedicine. Recently, there has been an increase in data regarding the involvement of not only dopaminergic neurons of the substantia nigra but also astrocytes in the pathogenesis of PD. Cell models based on induced pluripotent stem cells (iPSCs) and their differentiated derivatives are a useful tool for studying the contribution and interaction of these two cell types in PD. Here, we generated two iPSC lines, ICGi034-B and ICGi034-C, by reprogramming peripheral blood mononuclear cells of a patient with a heterozygous mutation c.1226A>G (p.N370S) in the GBA1 gene by non-integrating episomal vectors encoding OCT4, KLF4, L-MYC, SOX2, LIN28, and mp53DD. The iPSC lines demonstrate the expression of pluripotency markers and are capable of differentiating into three germ layers. We differentiated the ICGi034-B and ICGi034-C iPSC lines into astrocytes. This resulting cell model can be used to study the involvement of astrocytes in the pathogenesis of GBA-associated PD.

Oxidative stress monitoring in iPSC-derived motor neurons using genetically encoded biosensors of H2O2.

Oxidative stress plays an important role in the development of neurodegenerative diseases, being either the initiator or part of a pathological cascade that leads to the neuron's death. Genetically encoded biosensors of oxidative stress demonstrated their general functionality and overall safety in various systems. However, there is still insufficient data regarding their use in the research of disease-related phenotypes in relevant model systems, such as human cells. Here, we establish an approach for monitoring the redox state of live motor neurons with SOD1 mutations associated with amyotrophic lateral sclerosis. Using CRISPR/Cas9, we insert genetically encoded biosensors of cytoplasmic and mitochondrial H2O2 in the genome of induced pluripotent stem cell (iPSC) lines. We demonstrate that the biosensors remain functional in motor neurons derived from these iPSCs and reflect the differences in the stationary redox state of the neurons with different genotypes. Moreover, we show that the biosensors respond to alterations in motor neuron oxidation caused by either environmental changes or cellular stress. Thus, the obtained platform is suitable for cell-based research of neurodegenerative mechanisms.

Generation of induced pluripotent stem cell line, ICGi034-A, by reprogramming peripheral blood mononuclear cells from a patient with Parkinson's disease associated with GBA mutation.

Mutation in the glucocerebrosidase encoding gene (GBA) is one of the most frequent genetic cause of Parkinson's disease. ICGi034-A induced pluripotent stem cell (iPSC) line obtained by reprogramming peripheral blood mononuclear cells (PBMCs) of a patient with heterozygous c.1226A > G (p.N370S) mutation in the GBA gene can be used for studying the fundamental mechanisms of the pathogenesis of GBA-associated Parkinson's disease, and for potential drug screening. The iPSCs express pluripotency markers (NANOG, SSEA4, TRA-1-60, OCT4, SOX2), have a normal karyotype, and are capable of producing derivatives of three germ layers.

New Hybrid Compounds Combining Fragments of Usnic Acid and Thioether Are Inhibitors of Human Enzymes TDP1, TDP2 and PARP1.

Tyrosyl-DNA phosphodiesterase 1 (TDP1) catalyzes the cleavage of the phosphodiester bond between the tyrosine residue of topoisomerase 1 (TOP1) and the 3' phosphate of DNA in the single-strand break generated by TOP1. TDP1 promotes the cleavage of the stable DNA-TOP1 complexes with the TOP1 inhibitor topotecan, which is a clinically used anticancer drug. This article reports the synthesis and study of usnic acid thioether and sulfoxide derivatives that efficiently suppress TDP1 activity, with IC50 values in the 1.4-25.2 µM range. The structure of the heterocyclic substituent introduced into the dibenzofuran core affects the TDP1 inhibitory efficiency of the compounds. A five-membered heterocyclic fragment was shown to be most pharmacophoric among the others. Sulfoxide derivatives were less cytotoxic than their thioester analogs. We observed an uncompetitive type of inhibition for the four most effective inhibitors of TDP1. The anticancer effect of TOP1 inhibitors can be enhanced by the simultaneous inhibition of PARP1, TDP1, and TDP2. Some of the compounds inhibited not only TDP1 but also TDP2 and/or PARP1, but at significantly higher concentration ranges than TDP1. Leader compound 10a showed promising synergy on HeLa cells in conjunction with the TOP1 inhibitor topotecan.

Genetically Encoded Fluorescent Biosensors for Biomedical Applications.

One of the challenges of modern biology and medicine is to visualize biomolecules in their natural environment, in real-time and in a non-invasive fashion, so as to gain insight into their physiological behavior and highlight alterations in pathological settings, which will enable to devise appropriate therapeutic strategies. Genetically encoded fluorescent biosensors constitute a class of imaging agents that enable visualization of biological processes and events directly in situ, preserving the native biological context and providing detailed insight into their localization and dynamics in cells. Real-time monitoring of drug action in a specific cellular compartment, organ, or tissue type; the ability to screen at the single-cell resolution; and the elimination of false-positive results caused by low drug bioavailability that is not detected by in vitro testing methods are a few of the obvious benefits of using genetically encoded fluorescent biosensors in drug screening. This review summarizes results of the studies that have been conducted in the last years toward the fabrication of genetically encoded fluorescent biosensors for biomedical applications with a comprehensive discussion on the challenges, future trends, and potential inputs needed for improving them.

Mild phenotype of knockouts of the major apurinic/apyrimidinic endonuclease APEX1 in a non-cancer human cell line.

The major human apurinic/apyrimidinic (AP) site endonuclease, APEX1, is a central player in the base excision DNA repair (BER) pathway and has a role in the regulation of DNA binding by transcription factors. In vertebrates, APEX1 knockouts are embryonic lethal, and only a handful of knockout cell lines are known. To facilitate studies of multiple functions of this protein in human cells, we have used the CRISPR/Cas9 system to knock out the APEX1 gene in a widely used non-cancer hypotriploid HEK 293FT cell line. Two stable knockout lines were obtained, one carrying two single-base deletion alleles and one single-base insertion allele in exon 3, another homozygous in the single-base insertion allele. Both mutations cause a frameshift that leads to premature translation termination before the start of the protein's catalytic domain. Both cell lines totally lacked the APEX1 protein and AP site-cleaving activity, and showed significantly lower levels of the APEX1 transcript. The APEX1-null cells were unable to support BER on uracil- or AP site-containing substrates. Phenotypically, they showed a moderately increased sensitivity to methyl methanesulfonate (MMS; ~2-fold lower EC50 compared with wild-type cells), and their background level of natural AP sites detected by the aldehyde-reactive probe was elevated ~1.5-2-fold. However, the knockout lines retained a nearly wild-type sensitivity to oxidizing agents hydrogen peroxide and potassium bromate. Interestingly, despite the increased MMS cytotoxicity, we observed no additional increase in AP sites in knockout cells upon MMS treatment, which could indicate their conversion into more toxic products in the absence of repair. Overall, the relatively mild cell phenotype in the absence of APEX1-dependent BER suggests that mammalian cells possess mechanisms of tolerance or alternative repair of AP sites. The knockout derivatives of the extensively characterized HEK 293FT cell line may provide a valuable tool for studies of APEX1 in DNA repair and beyond.

The Cutting Edge of Disease Modeling: Synergy of Induced Pluripotent Stem Cell Technology and Genetically Encoded Biosensors.

The development of cell models of human diseases based on induced pluripotent stem cells (iPSCs) and a cell therapy approach based on differentiated iPSC derivatives has provided a powerful stimulus in modern biomedical research development. Moreover, it led to the creation of personalized regenerative medicine. Due to this, in the last decade, the pathological mechanisms of many monogenic diseases at the cell level have been revealed, and clinical trials of various cell products derived from iPSCs have begun. However, it is necessary to reach a qualitatively new level of research with cell models of diseases based on iPSCs for more efficient searching and testing of drugs. Biosensor technology has a great application prospect together with iPSCs. Biosensors enable researchers to monitor ions, molecules, enzyme activities, and channel conformation in live cells and use them in live imaging and drug screening. These probes facilitate the measurement of steady-state concentrations or activity levels and the observation and quantification of in vivo flux and kinetics. Real-time monitoring of drug action in a specific cellular compartment, organ, or tissue type; the ability to screen at the single-cell resolution; and the elimination of the false-positive results caused by low drug bioavailability that is not detected by in vitro testing methods are a few of the benefits of using biosensors in drug screening. Here, we discuss the possibilities of using biosensor technology in combination with cell models based on human iPSCs and gene editing systems. Furthermore, we focus on the current achievements and problems of using these methods.

New Hybrid Compounds Combining Fragments of Usnic Acid and Monoterpenoids for Effective Tyrosyl-DNA Phosphodiesterase 1 Inhibition.

Usnic acid (UA) is a secondary metabolite of lichens that exhibits a wide range of biological activities. Previously, we found that UA derivatives are effective inhibitors of tyrosyl-DNA phosphodiesterase 1 (TDP1). It can remove covalent complex DNA-topoisomerase 1 (TOP1) stabilized by the TOP1 inhibitor topotecan, neutralizing the effect of the drugs. TDP1 removes damage at the 3' end of DNA caused by other anticancer agents. Thus, TDP1 is a promising therapeutic target for the development of drug combinations with topotecan, as well as other drugs for cancer treatment. Ten new UA enamino derivatives with variation in the terpene fragment and substituent of the UA backbone were synthesized and tested as TDP1 inhibitors. Four compounds, 11a-d, had IC50 values in the 0.23-0.40 µM range. Molecular modelling showed that 11a-d, with relatively short aliphatic chains, fit to the important binding domains. The intrinsic cytotoxicity of 11a-d was tested on two human cell lines. The compounds had low cytotoxicity with CC50 ≥ 60 µM for both cell lines. 11a and 11c had high inhibition efficacy and low cytotoxicity, and they enhanced topotecan's cytotoxicity in cancerous HeLa cells but reduced it in the non-cancerous HEK293A cells. This "protective" effect from topotecan on non-cancerous cells requires further investigation.

Comparative Metabolomic Profiling of Rat Embryonic and Induced Pluripotent Stem Cells.

Metabolomic profiles of somatic cells, embryonic stem cells (ESCs), and induced pluripotent stem cells (iPSCs) reflect their metabolic phenotypes. The comparative study of metabolomes of these cells is important for understanding the differences in metabolism between somatic and pluripotent cells, and also the possible differences between ESCs and iPSCs. Here, we performed for the first time the metabolomic analysis of rat ESCs, iPSCs, and embryonic fibroblasts (EFs) at both quantitative and semi-quantitative levels using NMR spectroscopy and liquid chromatography with mass spectrometric detection, respectively. The total of 106 metabolites has been identified, and the concentrations of 51 compounds have been measured. It is found that the reprogramming of rat EFs into iPSCs affects virtually all metabolic pathways and causes drastic changes in the cell metabolomic profile. The difference between ESCs and iPSCs is much less pronounced: the concentrations of the majority of metabolites in ESCs and iPSCs are similar, and significant differences were observed for only several compounds, including adenosine, cysteic acid, glycerophosphoglycerol, inositol phosphate, glucose, myo-inositol, phosphoserine, xanthosine, guanosine. The observed differences between the metabolomic compositions of ESCs and iPSCs do not influence the pluripotent ability of iPSCs. Graphical Abstract.

Design, Synthesis, and Biological Investigation of Novel Classes of 3-Carene-Derived Potent Inhibitors of TDP1.

Two novel structural types of tyrosyl-DNA phosphodiesterase 1 (TDP1) inhibitors with hexahydroisobenzofuran 11 and 3-oxabicyclo [3.3.1]nonane 12 scaffolds were discovered. These monoterpene-derived compounds were synthesized through preliminary isomerization of (+)-3-carene to (+)-2-carene followed by reaction with heteroaromatic aldehydes. All the compounds inhibit the TDP1 enzyme at micro- and submicromolar levels, with the most potent compound having an IC50 value of 0.65 µM. TDP1 is an important DNA repair enzyme and a promising target for the development of new chemosensitizing agents. A panel of isogenic clones of the HEK293FT cell line knockout for the TDP1 gene was created using the CRISPR-Cas9 system. Cytotoxic effects of topotecan (Tpc) and non-cytotoxic compounds of the new structures were investigated separately and jointly in the TDP1 gene knockout cells. For two TDP1 inhibitors, 11h and 12k, a synergistic effect was observed with Tpc in the HEK293FT cells but was not found in TDP1 -/- cells. Thus, it is likely that the synergistic effect is caused by inhibition of TDP1. Synergy was also found for 11h in other cancer cell lines. Thus, sensitizing cancer cells using a non-cytotoxic drug can enhance the efficacy of currently used pharmaceuticals and, concomitantly, reduce toxic side effects.

Noncoding RNAs in the Regulation of Pluripotency and Reprogramming.

Pluripotent stem cells have great potential for developmental biology and regenerative medicine. Embryonic stem cells, which are obtained from blastocysts, and induced pluripotent stem cells, which are generated by the reprogramming of somatic cells, are two main types of pluripotent cells. It is important to understand the regulatory network that controls the pluripotency state and reprogramming process. Various types of noncoding RNAs (ncRNAs) have emerged as substantial components of regulatory networks. The most studied class of ncRNAs in the context of pluripotency and reprogramming is microRNAs (miRNAs). In addition to canonical microRNAs, other types of small RNAs with miRNA-like function are expressed in PSCs. Another class of ncRNAs, long ncRNAs, are also involved in pluripotency and reprogramming regulation. Thousands of ncRNAs have been annotated to date, and a significant number of the molecules do not have known function. In this review, we briefly summarized recent advances in this field and described existing genome-editing approaches to study ncRNA functions.

Impact of Xist RNA on chromatin modifications and transcriptional silencing maintenance at different stages of imprinted X chromosome inactivation in vole Microtus levis.

In vole Microtus levis, cells of preimplantation embryo and extraembryonic tissues undergo imprinted X chromosome inactivation (iXCI) which is triggered by a long non-coding nuclear RNA, Xist. At early stages of iXCI, chromatin of vole inactive X chromosome is enriched with the HP1 heterochromatin-specific protein, trimethylated H3K9 and H4K20 attributable to constitutive heterochromatin. In the study, using vole trophoblast stem (TS) cells as a model of iXCI, we further investigated chromatin of the inactive X chromosome of M. levis and tried to find out the role of Xist RNA. We demonstrated that chromatin of the inactive X chromosome in vole TS cells also contained the SETDB1 histone methyltransferase and KAP1 protein. In addition, we observed that Xist RNA did not contribute significantly to maintenance of X chromosome inactive state during iXCI in vole TS cells. Xist repression affected neither transcriptional silencing caused by iXCI nor maintenance of trimethylated H3K9 and H4K20 as well as HP1, KAP1, and SETDB1 on the inactive X chromosome. Moreover, the unique repertoire of chromatin modifications on the inactive X chromosome in vole TS cells could be disrupted by a chemical compound, DZNep, and then restored even in the absence of Xist RNA. However, Xist transcript was necessary for recruitment of an additional repressive histone modification, trimethylated H3K27, to the inactive X chromosome during vole TS cell differentiation.

Genome-wide profiling and differential expression of microRNA in rat pluripotent stem cells.

MicroRNAs (miRNAs) constitute a class of small noncoding RNAs that plays an important role in the post-transcriptional regulation of gene expression. Much evidence has demonstrated that miRNAs are involved in regulating the human and mouse pluripotency. Nevertheless, to our knowledge, miRNAs in the pluripotent stem cells of one of the most commonly used model organisms - the Rattus norvegicus have not been studied. In the present study, we performed deep sequencing of small RNA molecules in the embryonic fibroblasts, embryonic stem cells, and induced pluripotent stem cells of laboratory rats. Bioinformatics analysis revealed 674 known miRNAs and 394 novel miRNA candidates in all of the samples. Expression of known pluripotency-associated miRNAs, such as the miR-290-295 and miR-183-96-182 clusters as well as members of the miR-200 family, was detected in rat pluripotent stem cells. Analysis of the targets of differentially expressed known and novel miRNAs showed their involvement in the regulation of pluripotency and the reprogramming process in rats. Bioinformatics and systems biology approaches identified potential pathways that are regulated by these miRNAs. This study contributes to our understanding of miRNAs in the regulation of pluripotency and cell reprogramming in the laboratory rat.

Modern Genome Editing Technologies in Huntington's Disease Research.

The development of new revolutionary technologies for directed gene editing has made it possible to thoroughly model and study NgAgo human diseases at the cellular and molecular levels. Gene editing tools like ZFN, TALEN, CRISPR-based systems, NgAgo and SGN can introduce different modifications. In gene sequences and regulate gene expression in different types of cells including induced pluripotent stem cells (iPSCs). These tools can be successfully used for Huntington's disease (HD) modeling, for example, to generate isogenic cell lines bearing different numbers of CAG repeats or to correct the mutation causing the disease. This review presents common genome editing technologies and summarizes the progress made in using them in HD and other hereditary diseases. Furthermore, we will discuss prospects and limitations of genome editing in understanding HD pathology.

Dynamic properties of SOD1 mutants can predict survival time of patients carrying familial amyotrophic lateral sclerosis.

One of the reasons for the death of motor neurons of the brain and spinal cord in patients with amyotrophic lateral sclerosis is known to be formation of subcellular protein aggregates that are caused by mutations in the SOD1 gene. Patient survival time was earlier shown to have limiting correlation with thermostability change of SOD1 mutant forms of patients' carriers. We hypothesized that aggregation of mutant SOD1 may occur not only due to the protein destabilization, but through formation of novel interatomic bonds which stabilize "pathogenic" conformations of the mutant as well. To estimate these effects in the present paper, we performed statistical analysis of occupancy of intramolecular hydrogen bonds, hydrogen bonds between the protein and water molecules, and water bridges with use of molecular dynamics simulation for 38 mutant SOD1 forms. Multiple regression model based on these kinds of bonds demonstrated correlation with patient survival time significantly better (R = .9, p-value < 10-11) than the thermostability of SOD1 mutants only. It was shown that the occupancy of intramolecular hydrogen bonds between amino acid residues is a key determinant (R = .89, p-value < 10-10) in predicting patients' survival time.